Impurity resonance states in noncentrosymmetric superconductor CePt 3 Si: a probe 

for Cooper-pairing symmetry 
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Motivated by the recent discovery of noncentrosymmetric superconductors, such as CePtzSi, 
CeRhSii and CelrSiz, we investigate theoretically the impurity resonance states with coexisting 
s- and p-wave pairing symmetries. Due to the nodal structure of the gap function, we find single 
nonmagnetic impurity- induced resonances appearing in the local density of state (LDOS). In par- 
ticular, we analyze the evolution of the local density of states for coexisting isotropic s-wave and 
p-wave superconducting states and compare with that of anisotropic s-wave and p-wave symmetries 
of the superconducting gap. Our results show that the scanning tunneling microscopy can shed light 
on the particular structure of the superconducting gap in non-centrosymmetric superconductors. 
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INTRODUCTION 



Recent discoveries of superconductivity in the sys- 
tems that posses a lack of inversion symmetry such 
as CePt^Si^ with T c ~ 0.75K and more recently 
CeRhSi 3 2 , CelrSi^, Li(Pdi- x , Pt x ) a B±, UIt&, Y 2 C 3 6 
have raised an interest in the theoretical investigation 
of superconductivity in these systems. Among inter- 
esting questions the most important one concerns the 
underlying symmetry of the superconducting order pa- 
rameter. In particular, in all these materials, there is 
a nonzero potential gradient W averaged in the unit 
cell due to lack of inversion symmetry which results in 
the anisotropic spin-orbit interaction. Its general form 
can be determined by a group theoretical argument 7 
and, as it has been found, leads to many interest- 
ing propertie o 8 ' 9 ' 10 ! 11 ] 12 ! 13 ] 14 . For example, on general 
grounds there is a mixing of the spin-singlet and spin 
triplet superconducting states due to the lack of inver- 
sion. In CePtaSi the pairing symmetry has been stud- 
ied theoretically ^ 10 i 11 ' 12 i 13 and it is believed that the 
s + p-wave superconducting state is realized. Frigeri et 
aZjii pointed out that the spin-orbit interaction could 
determine the direction of the d- vector as d||i (I is the 
vector of the Rashba spin-orbit coupling) for which the 
highest transition temperature was obtained. A micro- 
scopic calculation with the detailed structure of the Fermi 
surface 1 — seems to confirm that the s+p wave state is the 
most probable one. The experimental studies of the tem- 
perature dependencies of the spin-lattice relaxation 15 , 
the magnetic penetration deptbiS, and the thermal con- 
ductivity measurements^ 7 " are also consistent with this 
conjecture. 

It is known that the non-magnetic as well as the mag- 
netic impurities in the conventional and unconventional 
superconductors already have been proven to be a useful 
tool to distinguish between various symmetries of the su- 
perconducting state 18 . For example, in the conventional 
isotropic s-wave superconductor the single magnetic im- 



purity induced resonance state is located at the gap edge, 
which is known as Yu-Shiba-Rusinov stated. In the 
case of unconventional superconductor with d x 2_ y 2-wa,ve 
symmetry of the superconducting state the non-magnetic 
impurity-induced bound state appears near the Fermi en- 
ergy as a hallmark of d x 2_ y 2-w&ve pairing symmetry 2 ^. 
The origin of this difference is understood as being due 
to the nodal structure of two kinds of SC order: in the 
d x 2_ y 2-wave case the phase of Cooper-pairing wave func- 
tion changes sign across the nodal line which yields finite 
density of states below the superconducting gap, while in 
the isotropic s-wave case the density of states is gapped 
up to energies of about Ao and thus the bound state can 
appear only at the gap edge. In principle the formation 
of the impurity resonance states can also occur in uncon- 
ventional superconductors if the nodal line or point does 
not exist at the Fermi surface of a superconductor like 
it occurs for isotropic nodeless p-wave and/or d x + id y - 
wave superconductors for the large value of the potential 
strength 21 . Therefore, STM measurements of the im- 
purity states can provide important messages about the 
pairing symmetry in the revelent systems. In the non- 
centrosymmetric superconductor with the possible coex- 
istence of s-wave and p-wave pairing symmetry, it is very 
interesting to see what is the nature of the impurity state, 
and whether a low energy resonance state can still oc- 
cur around the impurity through changing the dominant 
role played by each of the pairing components. Previ- 
ously the effect of the non-magnetic impurity scattering 
has been studied in the non-centrosymmetric supercon- 
ductors with respect to the suppression of T<~ and the 
behavior of the upper critical field 2 ^.. 

In this paper we investigate theoretically the impurity 
resonance states where both s-wave and p-wave Cooper- 
pairing coexist. Due to the nodal structure of gap func- 
tion as a result of the interference between the spin triplet 
and the spin singlet components of the superconduct- 
ing order parameters, we find that a single nonmagnetic 
impurity-induced resonance state appears in the local 
density of state. In particular, we analyze the evolu- 
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tion of the local density of states for coexisting isotropic 
s-wave and p-wave superconducting states and compare 
with that of anisotropic s-wave and p-wave symmetries 
of the superconducting gap. Our results show that the 
scanning tunneling microscopy can shed light on the 
particular structure of the superconducting gap in non- 
centrosymmetric superconductors. 



II. THE MODEL AND T-MATRIX 
FORMULATION 

Theoretical models of the superconducting state in 
CePtaSi are based upon the existence of a Rashba type 
spin-orbit coupling (RSOC)^. Therefore, following previ- 
ous consideration^ we start from an single orbital model 
with RSOC 



H = 



E 



; E 



gk • OSs'Cl Cks 



(1) 



where c^ s (cks) is the fermion creation (annihilation) op- 



erator with spin s and momentum k. 
tight-binding energy dispersion 



Here, £k is the 
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FIG. 1: (color online) The calculated Fermi surface using the 
Eq. ((T| and the spin-orbit coupling constant a = 0.3t. 



£k = 2t(cos(k x ) + cos(ky)) + 4ti cos(k x ) cos(ky) 

+ 2t 2 (cos(2k x ) + cos(2k y )) 

+ [2*3 + 4* 4 (cos(fe x ) + cos(fc y )) 

+ 4t 5 (cos(2fc x ) + cos(2k y ))] cos(k z ) 

+ 2* 6 cos(2fc z ) - /i (2) 

which reproduces the so-called /3-band of CePt^Si as ob- 
tained from the band structure calculation o 7 ' 13 . The elec- 
tron hopping parameters are (*, *i, *2, *3, *4, *5, *6, n) — 
(1,-0.15,-0.5,-0.3,-0.1,-0.09,-0.2,1.75) and the 
electron density per site n is used to determine the chem- 
ical potential^ 3 -. 

The second term of Eq.([T]) is the RSOC interac- 
tion where a denotes the coupling constant and the 
vector function gk is assumed in the following form 
gk = (— sink y ,smk x ,0). This term removes the usual 
Kramers degeneracy between the two spin states at a 
given k, and leads to a quasiparticle dispersion ek = 

£k ± a|gk| with |g k | = Jgkl +gky +Skl, splitting the 
Fermi surface (FS) into two sheets. Based on the above 
hopping parameters and RSOC constant a — 0.3*, the 
resulting FS is shown in FigfTJ where the main character- 
istic features of the FS has been successfully reproduced 7 . 

In the superconducting state, the presence of RSOC 
breaks the parity and, therefore, mixes the singlet (even 
parity) and triplet (odd parity) Cooper-pairing states. 
A full symmetry analysi ng 13 shows that s-wave pairing 
A s = Ao(cos(fc x ) + cos(fcy)) and a p-wave triplet pairing 
state with order parameter dk parallel to the gk vec- 
tor, dk = (iogk are able to coexist. Following previous 
estimations^ 3 - we have taken the odd parity component 



dk = do(— sin k y , sin k x , 0). Then the mean field BCS 
Hamiltonian for this system has the matrix form 
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idky). The inverse of the single-particle Green's function 
is defined as 



g 1 (k, iu n ) = iuj n I - i? k , 



(4) 



where / is the 4x4 identity matrix. Taking the inverse 
of Eq. © we find 



g(k,iuj n ) = 
where 

G{k,iw n ) 



G(k, iuj n ) F(k,iu„) 
Ft(k,iw n ) -G*(-k,-iw„) 



(5) 



E 

r=±l 

F(k,ico n ) = Y, 



1 + r(g k • a-) 



G T (k,ioj n ), (6) 



1 + 7-(g k • a) . 



r=±l 



and 



G T (k, iwn) 
F T (k, iw n ) 



A r 



[ILO,,-, 



ia y F T (k,iuj n ), (7) 
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Here, the single-particle excitation energy is 
with 



£k + ralgkl; A T 



rldL 



(10) 



(11) 



and the unit vector is gk = gk/|gk|- 

For completion the equations above have to be supple- 
mented by the self-consistency equation that determines 
the symmetry of the superconducting gap and the super- 
conducting transition temperature. For the sake of sim- 
plicity and also because this is not critical for our further 
analysis we consider the superconducting order param- 
eter as a given parameter. At the same time, recent 
studies based on the helical spin fluctuation mediated 
Cooper-pairing find two stable superconducting phases 
with either dominantly s +p-wave or p + d + /-wave sym- 
metry of superconducting order paramete r 13 i 24 ' 25 . In the 
following we adopt the former one for our calculation. 

The next step is to obtain the Green's function in the 
presence of a single impurity site. The impurity scatter- 
ing is given by 



(12) 



where without loss of generality we have taken a single- 
site nonmagnetic impurity of strength Uq located at 
the origin, = 0. Then the site dependent Green's 
function can be written in terms of the T-matrix 
formulatio n 21 ' 26 i 27 as 

C(i,j;iu} n ) = (o(i-j;iUn) 

+ Co(i,iw n )T(iuj n Xo(j,icOn), (13) 



where 



U p 3 



(14) 



1 - U p3(o(0, 0;iw„) 
( {iJ;iuj n ) = ^^e lk - Ri J 5 (fc,iw„), (15) 

k 

with pi being the Pauli spin operator, and Ri is the lattice 
vector, Ry = Ri Rj. Finally, the local density of state 
which can be measured in the STM experiment has been 
obtained as 

N(r, uj) = —— Im(ii(r, r;u; + irj), (16) 

i 

where rj denotes an infinitely small positive number. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. The density of state 

Before considering the effect of the impurity it is useful 
to analyze first the density of state (DOS) in the super- 



conducting state, which is expressed as, 

p(u) = ImVj„(k,w) (17) 

7T » — ' 

i,k 

As we already have mentioned above it is not neces- 
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FIG. 2: (color online) The evolution of the local density of 
states for various ratio between coexisting isotropic s-wave 
and p-wave Cooper-pairing state. The left and right panels 
refer to the different values of the damping constant V. The 
dashed and the dotted curve denote the contribution of the 
different bands and the straight curve refers to the total den- 
sity of states. The parameters of the gaps and the damping 
r are given in terms of hopping integral t. 

sary to calculate the magnitude of the gap functions self- 
consistently since we are mainly interested in the qual- 
itative properties arising from the gap structure. We 
first consider the situation when the s-wave part of the 
total superconducting gap is momentum independent, 
A s = Aq. In FigfJ] we show the evolution of the den- 
sity of state for positive frequencies for various values of 
the s-wave component of the superconducting gap. In 
particular, for zero value of the s-wave component the 
superconducting gap is purely determined by the p-wave 
superconducting gap with the point node at the Fermi 
surfaces of the corresponding bands at (k x = 0, k y = 0). 
This gap structure is the same for both bands splitted 
by the spin-orbit coupling. With increasing value of the 
isotropic s-wave gap one finds that the total supercon- 
ducting gap in one of the bands increasing with the to- 
tal superconducting gap A s + |dk| while it decreases ef- 
fectively for the other band for which the total gap is 
A s — |dk|- Once both s— wave and p-wave superconduct- 
ing gaps are the same, the accidental node forms at one 
of the band and the behavior of the density of states 
changes to a linear at low energy reflecting the formation 
of the line of node. We further note that density of states 
shows only slight electron-hole asymmetry. 
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In Fig. [3] we show a similar evolution of the den- 
sity of states, however, now the s-wave component of 
the superconducting gap is momentum dependent, A s = 
Aofcos/cz + cos/cy) = Ao7i<r^. Interestingly enough, here 
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FIG. 4: (color online) The LDOS for coexisting isotropic s- 
wave and p-wave Cooper pairing states for various ratio of the 
parameters. The straight (red) curves refer to the calculated 
density of states without impurity and the dashed (green) 
curves refer to the LDOS at the (0, 1,0) position. Here, we 
use Uo = 5t. 



FIG. 3: (color online) The evolution of the local density 
of states for various ratio between coexisting anisotropic s- 
wave and p-wave Cooper-pairing state. The left and right 
panels refer to the different values of the damping constant 
r. The dashed and the dotted curve denote the contribution 
of the different bands and the straight curve refers to the 
total density of states. The parameters of the gaps and the 
damping F are given in terms of hopping integral t. 

the node in the density of states forms already when 
the p-wave superconducting gap component is zero (see 
also Fig(5j) and is the result of the initial momentum 
structure of the s-wave superconducting gap that yields 
point nodes on the Fermi surface. This is unique to the 
anisotropic s-wave superconducting gap. By introducing 
the interference between s-wave and p-wave gap the po- 
sition of the node is shifted to the different points of the 
Brillouin Zone, however, here the nodal structure of the 
superconducting gap is not a result of the interference 
effect between p-wave and s-wave of the superconduct- 
ing gap but arises already in the pure anisotropic s-wave 
symmetry and shifted by introducing the moderate com- 
ponent of the p-wave gap. 




FIG. 5: (color online) The LDOS for coexisting anisotropic 
s-wave (A)and p-wave Cooper pairing states for various ratio 
of the parameters. The straight (red) curves refer to the cal- 
culated density of states without impurity and the dashed 
(green) curves refer to the LDOS at the (0,1,0) position. 
Here, we use Uo — 5t. 



B. Impurity resonance states 

In view of complicated band structure arising in 
CePtaSi from the Rashba spin-orbit coupling and the 
corresponding interference effect for the superconducting 
gap the density of states in a clean case that can be ac- 
cessed by the tunneling experiments cannot give a precise 
information on the exact structure of the superconduct- 
ing gap in the non-centrosymmetric superconductors. At 



the same time, an introduction of the non-magnetic im- 
purity can give additional important information on the 
symmetry of the superconducting gap in such a material. 
In terms of Eq.(16), the T-matrix can be written as 

T~ 1 (iu n ) = IV 1 -p3Co(0,0;icj n ), (18) 

and the position of the impurity resonant state is given by 
det T -1 = 0. We first study the situation of the isotropic 
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s-wave superconducting gap coexisting with p-wave. In 
FigH] we show the calculated density of states without 
impurity and also the local density of states with an im- 
purity on the nearest neighbor site (0, 1,0). Without the 
s-wave component the density of states shows the forma- 
tion of the impurity induced resonant bound states that 
appear symmetrically in energy at the positive and neg- 
ative sides of the LDOS. Clearly these resonant bound 
states arise due to unconventional nature of the p-wave 
superconducting gap and the nodal points at the Fermi 
surface. One clearly sees that upon increasing of the iso- 
topic s-wave contribution the bound state shifts towards 
the edge of the superconducting gap implying the zero 
density of states for energies lower than Ao . 

In FiglU we show the corresponding local density of 
states for the coexisting anisotropic s-wave and p-wave 
superconducting gaps. In the present case, for any value 
of the s-wave and p-wave gap there are nodal points at 
the Fermi surface resulting either from the internal struc- 
ture of the anisotropic s-wave gap, point nodes from the 
p-wave state, or a nodal line at one of the bands that 
arises due to interference of the p-wave and s-wave gap. 
Therefore, the impurity induced bound state occurs for 



all ratios between the p-wave and s-wave gap. Note, 
that in case of pure anisotropic s-wave gap due to the 
nodal structure on both of the bands the impurity in- 
duced bound state becomes visible only for a very large 
values of the potential scattering strength Uq. 



IV. SUMMARY 

In summary, we have investigated theoretically the 
non-magnetic impurity induced resonance bound states 
in the superconductors without inversion symmetry us- 
ing as an example CePt^Si, which is believed to have 
a line node in the energy gap arising from the coexis- 
tence of s-wave and p-wave pairing symmetry. Analyz- 
ing the local density of states we find that the nodal 
structure of gap function, we find that a single nonmag- 
netic impurity-induced resonance states is highly prob- 
able in non-centrosymmetric superconductors. We show 
that further STM experiments may reveal the exact sym- 
metry of the superconducting gap in these systems. 

We thank Jun Chang for useful discussions. 
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